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ABSTRACT 


The  new  approach  to  structural  reliability  analysis  based  on 
order  statistics  and  the  expected  time  to  first  failure  in  a  fleet 
of  specified  magnitude,  that  has  been  first  proposed  in  a  previous 
report  AFML-TR-66-37  is  developed  in  more  detail  and  applied  to 
structures  subject  to  progressive  fatigue  damage.  It  is  shown  that 
the  concept  of  fatigue  sensitivity  developed  in  WADD  Tech .  Rep .  61- 
53  is  improved  by  relating  it  to  the  expected  time  to  first  failure 
rather  than  to  the  expected  time  to  failure.  The  modified  fatigue 
sensitivity  factor  becomes  an  effective  parameter  for  the  correla¬ 
tion  of  ultimate  load  and  fatigue  design  and  for  the  classification 
of  fatigue  sensitive  structures. 


This  abstract  is  subject  to  special  export  controls  and  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  the  Metals  and  Ceramics  Division  (MAM) ,  Air  Force  Materials 
Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio  45433. 
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I.  Introduction 


A  fundamental  difference  in  the  approach  to  the  life  and 
reliability  estimate  of  easily  maintainable,  multi-element  elec¬ 
tronic  systems  and  of  major  structures  the  performance  of  which 
depends  on  a  relatively  small  number  of  critical,  and  essential¬ 
ly  "non-maintainable"  elements  results  from: 

(a)  the  different  implication  of  the  concept  of  "failure", 

(b)  the  large  difference  in  the  possible  number  of  test 
replications  to  provide  an  estimate  of  "mean  life",  and 

(c)  the  basic  difference  in  the  procedures  of  reliability 
demonstration . 

In  electronic  systems  failure  in  operation  is  considered  a 
contingency,  reducible  by  a  program  of  periodic  replacement  of 
components  to  an  acceptable  condition  defined  by  a  small  number 
of  failures.  The  mean  time  to  failure  of  elements,  on  the  reli¬ 
able  estimate  of  which  the  effective  design  of  the  replacement 
program  depends,  is  obtained  from  tests  with  a  sufficiently  large 
number  of  replications.  For  a  major  structure  failure  is  not  a  con¬ 
tingency  but  an  emergency.  Since  the  consequences  of  even  a  single 
failure  are  usually  severe,  the  concept  of  "mean  life",  with  its 
tacit  implication  of  a  high  percentage  of  failed  structures  has 
no  relevance  to  design.  It  can  also  not  be  reliably  estimated 
from  the  single  test  or  the  very  small  number  of  test  replications 
that  are  technically  and  economically  feasible.  Therefore  the 
approach  to  reliability  demonstration  is  quite  different:  in 
electronic  systems  reliability  demonstration  of  the  whole  system  by 
computation  based  on  experimentally  demonstrated  component  reliability 
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is  acceptable,  while  even  the  crudest  estimate  of  the  reliability 
of  a  full-scale  structure  requires  a  life  test  of  this  structure 
itself,  supplemented  by  knowledge,  accumulated  in  multiple  tests 
of  large  structural  elements,  of  the  expected  range  of  scatter. 

This  difference  in  approach  must  be  reflected  in  a  difference 
in  methodology.  While  in  reliability  analysis  of  electronic 
systems  linear  functional  models  of  utmost  simplicity,  such  as 
the  simple  chain  or  parallel  chains,  form  the  substructure  upon 
which  probabilistic  models  of  different  levels  of  sophistication 
are  superimposed,  functional  models  representing  structural  failure 
are  by  the  nature  of  the  physical  process  highly  nonlinear  and  re¬ 
dundant.  Sophistication  in  the  superimposed  probability  model  would 
therefore  compound  the  difficulties  of  analysis  to  an  extent  that 
does  not  seem  warranted  by  the  obtainable  results  in  view  of  the 
severely  limited  number  of  experimental  data  available  for  their 
verification.  Since  under  such  conditions  reliability  demonstration 
must  rely  on  a  combination  of  full  scale  testing  and  performance 
during  early  life,  preferably  under  operational  conditions  of  more 
than  normal  severity,  the  expected  time  to  the  first  failure'*'  seems 
to  be  a  concept  that  is  more  relevant  to  the  life  estimate  and  the 
reliability  demonstration  of  major  structures,  in  the  design  of 
which  the  avoidance  of  even  a  single  catastrophic  failure  represents 
the  overriding  design  criterion,  than  the  conventional  concepts  of 
expected  (mean,  median  or  modal)  life  and  scatter  factor.  Obviously, 
failure  need  not  necessarily  be  defined  as  catastrophic,  but  can 
refer  to  a  designated  level  of  damage.  In  this  case  the  expected 
time  to  the  second  failure  as  well  as  the  expected  interval  between 
first  and  second  failures  become  useful  as  additional  characteristics 
from  which  not  only  the  "expected  life"  might  be  estimated,  but  which 
might  serve  to  verify  the  validity  of  various  possible  assumptions 
concerning  the  form  of  the  distribution  function  of  fatigue  life. 
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The  concepts  of  the  expected  time  to  the  first  failure  as  well  as 
that  of  the  interval  between  the  first  and  second  failures  or  that 
of  the  time  to  the  second  failure  also  have  the  practical  advantages 
that 

(i)  they  can  be  actually  demonstrated  in  operation,  which  is 
impossible  for  the  mean  time  to  failure, 

(ii)  they  are  strongly  affected  by  differences  between  different 
distributions  of  life  so  that  they  automatically  reflect  the  differences 
between  different  failure  mechanisms  which  lead  to  different  distrib¬ 
utions  of  times  to  failure,  and 

(iii) they  depend  on  the  size  of  the  group  of  structures  to  be 
considered . 

It  seems,  therefore,  to  be  of  practical  interest  to  study  the 
relation  between  these  concepts  and  the  conventional  reliability 
concepts  of  mean  time  to  failure  and  scatter  range,  because  of  the 
advantage  of  a  structural  reliability  analysis  based  on  order  statistics 
rather  than  on  conventional  statistics. 

II.  Distribution  Functions  of  Physical  Significance 

In  both  materials  tests  and  life  tests  of  structures  the  number 
of  test  replications  that  is  technically  and  economically  feasible 
precludes  the  derivation  of  the  distribution  function  of  fatigue 
lives  from  observations  by  methods  of  statistical  inference.  It  is 
therefore  necessary  to  devise,  on  the  basis  of  relevant  physical 
arguments,  probabilistic  models  which  produce  distribution  functions 
that  are  germane  to  the  physical  phenomenon,  so  that  their  use  can 
be  justified  in  the  extrapolation  from  small  numbers  of  test  results 
to  probability  levels  as  far  removed  from  the  range  of  test  results 
as  the  levels  that  are  significant  in  reliability  analysis.  However, 
it  should  be  realized  that  only  very  simple  physical  arguments  can 
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£>s  translated  into  probability  statements  which  result  in  a  charac¬ 
teristic  distribution  function. 

Thus,  for  instance,  the  argument,  based  on  observations  of  step¬ 
like  crack  propagation  under  cyclic  loading,  that  each  new  crack  ele¬ 
ment  can  be  considered  a  partial  static  rupture,  independent  of  pre¬ 
ceding  ruptures,  leads  to  the  conclusion  that  cracking  is  a  Poisson 
2 

process  .  The  probability  of  failure  associated  with  the  exceedance 
of  a  finite  number  r.  of  elementary  cracks  is  thus  given  as  the  sum  of 
a  finite  number  of  Poisson  terms 

r  -L/T  r 

P{x  s  r}  =  E  p(x)  =  «  E  ^r-  (t/T)X  =  P(t)  (2.1) 

o  o 


where  T  is  the  return  period  of  the  elementary  crack  process  of 
probability  p  =  T  ,  and  x  the  number  of  occurrences  in  a  given  time 
^•nbcrval  b .  The  time  to  failure  as  a  function  of  r  (mortality  function) 
is  obtained  by  differentiation  with  respect  to  t  of  P*=[l-P(t)] 

dP*/dt  =  (t/T) r  e  "t/T  (2.2) 

which,  with  x  ■  t,  a  ■  r  +  1,  X  =  T_1  and  r(a)  =  (o-l)J=r:  is  the 
Gamma  function 


p(x^  “  TToT  x(a  e  XX  (2.3) 

Alternatively,  if  it  is  assumed  that  an  elementary  crack  length 
e  is  exponentially  distributed,  the  momentary  crack  length  l  =  £Ei 
follows  a  Gamma  distribution^  for  the  crack  length  x  =  t  with  ^ 
a  “  k  and  X  =  1/e,  where  e  is  the  mean  elementary  crack  length, 
since  this  is  the  distribution  of  the  sum  of  k  random  variables, 
each  of  exponential  distribution.  Since  time  to  failure  can  be 
associated  with  the  time  to  attain  a  certain  crack  length,  since 
this  is  the  distribution  of  the  sum  of  k  random  variables,  each  of 
exponential  distribution.  Since  time  to  failure  can  be  associated 
with  the  time  to  attain  a  certain  crack  length,  both  physical  argu¬ 
ments  lead  to  Gamma  distributions  for  the  time  to  failure. 
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The  physical  assumption  of  an  increase  of  the  rate  of  cracking 
with  increasing  number  N  of  load  applications  and  thus  with  increasing 
time  t  can  be  translated  into  the  probabilistic  requirement  of  a  risk 
function  or  "hazard  rate",  increasing  with  increasing  t.  The  simplest 
form  of  this  function  satisfying  such  assumption 


a-1 


h  (t )  =  a  (-) 


a  >  1 


(2.4) 


where  v  is  a  measure  of  the  central  tendency  of  t  leads  to  the  extreme 
value  distribution 


P*  (x) 


- (x/v)a 

—  -6 


(2.5) 


with  x  =  t;  the  scale  parameter  a  is  an  inverse  measure  of  the  scatter 

of  x  or  t,  which  implies  that  the  more  rapidly  the  increase  of  the 

risk  function  the  narrower  the  associated  scatter  of  lives.  Eq.(2.5) 

represents  the  Third  asymptotic  distribution  function  of  the  shortest 
4 

lives  in  large  samples  of  well-behaved  statistical  populations  of 
unknown  distribution,  also  known  as  the  Weibull  distribution. 

Since  experiments  have  shown  that  a  certain  "incubation  period" 
t  ,  the  minimum  life  ,  is  needed  to  start  the  process  of  crack  propa¬ 
gation,  a  more  appropriate  formulation  of  the  hazard  rate 


h  (N)  =  a 


t-t 


a-1 


w-toJ 


a  >  1 


(2.6) 


leads  to  an  extreme  value  distribution  with  lower  limit 


x-x  ° 
o 


v-x 


P* (x)  =  e 


(2.7) 
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with  x  =  t  and  x  =  t  . 

o  o 

The  wide  use  of  the  Logarithmic  Normal  distribution  for  fatigue 
lives  is,  in  general,  based  on  arguments  of  expediency  (use  of  tables 
and  methods  of  estimation  developed  for  Normal  distribution) ,  rather 
than  physical  or  probabilistic  reasoning,  although  an  argument  that 
the  damaging  effect  of  each  load  application  is  proportional  to  the 
total  damage  produced  by  the  preceding  load  cycles  can  be  used  to 
justify  the  logarithmic  normal  distribution  of  fatigue  lives6.  The 
application  of  the  central  limit  theorem  to  a  product  of  statistical 
variables  provides  further  justification,  although,  physically,  the 
above  arguments  do  not  carry  conviction. 

The  hazard  functions  of  the  three  distribution  functions  are 
quite  dissimilar:  for  the  Gamma  distribution  it  increases  asymptotically 
toward  a  constant  value  associated  with  chance  failures,  for  the  Weibull 
distribution  it  increases  monotonicaly,  for  the  logarithmic  normal  it 
first  increases  at  decreasing  rate,  but  later  decreases  very  slowly 
towards  zero.  However,  an  experimental  demonstration  of  the  tendencies 
of  these  functions  is  just  as  impractical  as  a  demonstration  of  the 
forms  of  the  associated  reliability  or  frequency  functions  themselves 
(Fig.  1)  since  the  significant  differences  in  their  trend  only  appear 
in  the  central  range  of  lives. 

Certain  physical  assumptions  arising  from  considerations  of 
structural  failure  phenomena  lead  to  composite  distribution  functions 
of  various  types.  Thus,  for  instance,  the  assumption  that  fatigue 
failure  of  a  structure  subject  to  variable  load  intensities  is 
essentially  an  ultimate  load  (chance)  failure  of  the  structure  that 
has  been  progressively  weakened  by  propagating  fatigue  cracks  may 
be  expressed  by  replacing  the  constant  "return  period"  of  ultimate 
load  failure  v  =  1/Pf  in  the  exponential  distribution  characteristic 
for  chance  failures 
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P*  (x) 


(2.8) 


- (x/v) 

=  -e 

by  a  decreasing  function  v(x)  of  x  reflecting  the  progressive 
deterioration  by  fatigue  of  the  initial  ultimate  load  resistance  R. 

The  probability  of  ultimate  load  failure  of  the  undamaged  structure 
is  obtained  from  the  distribution  function  P ( v)  of  the  safety  factor 
v  =  R/S,  where  R  and  S  denote  respectively  the  resistance  and  the 

g 

load  considered  as  statistical  variables  ,  at  v  =  1>  or  P(v)=  P(1)=P  . 
Various  assumptions  could  be  made  concerning  the  form  of  the  function 
v(x).  Of  these  the  function 

v(x)  =  v  [l-(-)n]  >  o  with  n  >  1  (2.9) 

o  u 

seems  best  suited  to  reproduce  actual  conditions.  The  parameter  _u 
denotes  the  time  x  at  which  the  resistance  R  has  been  completely 
destroyed  by  fatigue  cracking.  Obviously  the  time  u  is  a  function 
of  the  operational  load  spectrum  that  causes  fatigue,  of  the  fatigue 
resistance  of  the  structure  which,  in  turn,  depends  on  the  fatigue 
performance  of  the  material  and  on  details  of  the  structural  design, 
and  on  the  maintenance,  repair  or  replacement  program  as  a  result  of 
which  the  time  u  is  periodically  extended.  Moreover  u  is  not  a  constant 
but  a  statistical  variable,  since  it  represents  the  fatigue  life  of 
the  structure  subject  to  the  normal  operational  load  spectrum  alone, 
independently  of  (extreme)  load  spectra  that  are  likely  to  produce 
actual  failure  (thunderstorm  gusts,  low-level  operations,  etc.).  In 
Eq.  (2.9)  this  variable  is  replaced  by  its  expected  (mean)  value.  The 
power  n  >  1  reproduces  the  observation  that  the  effect  of  fatigue 
deterioration  of  the  static  resistance  proceeds  at  first  very  slowly 
but  accelerates  sharply  as  x  approaches  u.  Introducing  Eq.  (2.9)  into 
(2.8)  the  distribution  function  is  obtained 


The  risk  function  associated  with  Eq.  (2.10)  follows  from 


H  H  (X/Vn} 

.  ,  d  .  .  d  o 

h(x)  =  -  inP*  =  -  - 

'  dx  dx 


1+  (n-1)  (^)n 

Cl—  (~ ) n ] — 2 


(2.11) 


Developing  the  exponent  of  Eq.  (2.10)  into  a  power  series 


n-0 


n  -1 
]  = 


(jj>  ^  + 


n 

<3>  + 


<  1 


(2.12) 


which  for  values  n  >  1  converges  rapidly  in  the  region  of  practical 
significance  (^)  <0.5,  Eq.  (2.10)  can  be  written  in  the  form 


P*(x) 


-  i<£> 

Cp  u 


n+1 


-  i(5> 

cp  U 


2n+l 


=  tt  -e 
m 


1 


mn+1 


(-) 

u 


TT  P*  (x) 
m  m 


(2.13) 


for  (x/u)  <  1  and  m  =  0,1,2,...,  where  the  "fatigue  ratio"  cp  =  (VQ/U) 

relates  the  return  period  of  ultimate  load  failure  to  the  expected 

time  to  complete  loss  of  resistance  caused  by  fatigue. 

The  function  P* (x)  in  Eq.  (2.13)  is  a  product  of  extreme  value 

distributions  P* (x)  of  the  type  of  Eq.  (2.5)  with  a  =  mn+1.  Since 
m  m 

a  is  inversely  proportional  to  the  standard  deviation  a  (logint)  =  t 
m  m  xu 


m 
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(2.14) 


6  =  rr/2.303a 

m  m 


it  follows  that  the  dispersion  of  the  functions  P  (x)  decreases  with 

in 

increasing  m. 

The  risk  function  associated  with  Eq.  (2.13) 


h(x) 


d 

dx 


inP* (x) 


.  mn+1 
m  u 


i  in 

I  (mn+1) (-) 
v  m  u 


(2.15) 


for  (G}<  is  an  increasing  function  of  x,  resulting  from  linear 
superposition  of  processes  of  progressive  damage.  For  large  values 
of  n  the  series  converges  fast  enough  to  justify  the  approximation 


h  (x)  =  i  +  (n+1) 
o 


(n+1) 


(2.16) 


which  represents  the  simple  superposition  of  a  chance  risk  and  a 
progressive  damage  risk  function  with  expected  life  =  u  ^"Vcp  >  u. 
P*  (x)  is  expressed  by  the  first  two  terms  of  Eq.  (2.13): 


-  f(x/u) 
P*  (x)  =  e  'P 


-  ~  (x/u) n+1 

9 


(2.17) 


Distribution  functions  of  the  general  type  of  Eq.  (2.13) 


P*  (x) 


* 


TT  P 


m  m 


(x) 


(2.18) 


arise  dn  the  reliability  analysis  of  multicomponent  structures  or  of 
structures  subject  to  a  multitude  of  independent  effects.  The  number 
m  of  product  terms  may  represent  the  number  of  components  surviving 
at  time  x  or  the  number  of  effects  acting  at  time  x  on  the  surviving 
structure . 
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HI-  Application  of  Order  Statistics 


Certain  order  statistics  ofseveral  distribution  functions  are 
evaluated  and  compared  with  the  expectation  that  the  results  will 
provide  a  basis  for  comparison  of  tests  and  ovservations  on  structures 
and  structural  parts  which  is  more  effective  than  conventional 
statistical  methods,  as  well  as  a  basis  for  the  establishment  of  a 
rational  and  administrable  procedure  of  reliability  demonstration. 

It  is  instructive  to  evaluate  first  the  order  statistics  for 


the  exponential  distribution  used  to  represent  chance  failures. 

The  expected  value  of  the  m-th  order  statistics  of  a  sample 
of  size  n  from  an  exponential  population  specified  by  Eq.  (2.8) 


and 


P  (x)  =  1  -  P*  (x)  =  1  -  .e 


/  \  1  -  (x/v) 

P(x)  =  -  -6 


-  (x/v) 


(3.1) 

(3.2) 


where  v  is  the  expected  mean  time  to  failure 


Hence 


and 


E  (x  ) 
m  n 


m 

v  E 
i=l 


1 

n-i+1 


(3.3) 


E(X1)  =  n  V;  E(X2)  = 


r  1  1  i 

[  -  +  - r  ]  V 

n  n-1 


(3.4) 


E(x2-x1)  -E(x2)  -  E(x1)  -sij  V  ,3  5) 

When  the  sample  size  is  moderate  or  large,  the  expected  times  to  first 
and  second  failure  are  quite  short,  much  shorter,  as  will  be  seen,  than 
these  times  for  any  of  the  considered  distributions  of  fatigue  life. 
This  is  an  important  fact  when  design  criteria  for  ultimate  load 
failure  are  compared  with  design  criteria  for  fatigue.  If  avoidance 
of  even  a  single  catastrophic  failure  is  a  design  requirement,  expected 
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-3 


times  to  first  failure  must  replace  the  mean  times  to  failure  as 
significant  design  criteria.  In  the  case  of  ultimate  load  failure 
the  requirement  of  an  expected  time  to  first  failure  of  x 1  =  100,000  hrs, 
in  a  fleet  of  n  =  1000  units  would  imply  a  design  mean  time  to  failure 

g 

for  the  fleet  of  v  =  10  hrs.,  producing  a  reliability  function 

g 

P*(x)  =  exp (-x/10  )  with  a  probability  of  failure  of  P  =  1  -  e-10 
-3  5  F 

10  at  x  =  x  =  10  .  The  expected  interval  between  the  first  two 

.  5 

consecutive  failures  is  about  10  hrs. 

For  a  distribution  function  P(x)  and  frequency  function  p(x)  the 
probability  that  no  failure  will  occur  in  the  interval  between  o  and  x 
and  one  of  a  sample  of  size  n  will  fail  in  the  interval  between  x  and 
(x+Ax)  is  given  by  the  expression 


p  (X^) -n[l-P  (x)  ]n  1  p  (x)  Ax 


(3.6) 


The  expected  time  to  the  first  failure  is  the  expectation  of  Eq.  (3.6), 
or 

EX, 


1  x[l-P(x)]n  ^(xjdx 


{i  -  n  1 

—  00 

Similarly,  the  time  to  the  second  failure 


(3.7) 


+00 


EX, 


=  n(n-l)  J  xP(x)  [l-P(x)  ]n  2p(x)dx  (3.8) 


The  expected  value  of  the  interval  between  the  first  and  second  failures 

(3.9) 


E(X2-X1)  =  EX2  -  EX 


The  variance  of  X,  is  obtained  from 


Var  X,  =  EX12  -  (EX^2 


where 


EX  =  n 

J. 

—  00 


x  [1-P(x)]n  1p(x)dx 


(3.10) 


(3.11) 
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(3.12) 


The  variance  of  X 2 

Var  X2 


(ex2)2 


where 

~f*°° 

EX  2  =  n  (n-1)  I*  x2P(x)ri-P(x)f  p(x)dx 

2  (3.13) 

The  coefficient  of  variation  of  P(x) 

V  =  /VS T  X/EX  (3.14) 

x 

9 

For  the  Third  Asymptotic  Extremal  Distribution  according  to  Eq. 

(2.5)  with  q 

p(x)  =  1  -  P*(x)  =  1  -  -e  (X/V)  (3.15) 

and 

p  (x )  =  (a/v)  (x/v)  «.  (3.16) 

evaluation  of  Eqs.  (3.7)  to  (3.13)  produces  the  expressions: 


EX^  =  vn  (1+1/a)  (3.17) 

EX12=  v2n"2/ar (1+2/a)  (3  18) 

EX2  =  vr  (1+1/a)  [n  (n-1)  1//a-(n-l)n  ^  ]  (3.19) 

7  2  —2/a  —2/ql-. 

EX2  =  v  T  (1+2/a)  [n  (n-1)  -(n-l)n  ]  (3.20) 


while 


and 


EX  =  vT (1+1/a) 

2  -.2 

Var  X=  [r  (1+2/a)  -  T  (1+1/a)]  v 


(3.21) 

(3.22) 
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Hence 


and  therefore 


-1/a 

EX1  =  n  EX 


2  -2/a  2 

EX1  =  n  EX 


-2/a 

Var  X^=  n  Var  X 


(3.17a) 


(3.18a) 


(3.23) 


o(X1)  =  n"1/ac(X)  (3.24) 

The  coefficient  of  variation 

vx  =  a(X1)/EX1  =  o  (X)  /EX  (3.25) 

is  therefore  identical  with  that  of  the  parent  population. 

Eqs.  (3.17a)  to  (3.24)  can  also  be  used  to  predict  the  parameters 
of  the  parent  population  from  observations  of  the  times  to  first 
failure  in  samples  of  size  n 


and 


EX  = 


1/a 

=  n  EX, 


a(X) 


1/a 

n 


o(Xx) 


(3.17b) 


(3.24a) 


If  Eq (2.7)  is  considered  with  the  lower  bound  (minimum  life) 

specified  as  a  fraction  of  the  characteristic  life  xq=  pv  so  that 

„  a 
rx-pv  -1 

Lv(i-p) J 

p*(x)  =  -6  (2.7a) 


the  transformation 

X  =  Y  +  pv 


(3.26) 
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reestablishes  the  form  of  Eq.  (3.15)  with  respect  to  the  variable  Y 


P(y) 


1-P* (y) 


- (y/w)a 

i— -e 


(3.27) 


where  w  =  v(l-p).  Hence  Eqs.  (3.17)  to  (3.20)  are  valid  with  respect 
to  the  variable  Y  and  the  following  expressions  can  be  written  for  the 
order  statistics  of  Eq.  (3.27): 

EY1  =  wn  1//ar  (1+1/a) 

or,  since 

\  =vpv 

the  expectation 


(3.28) 

(3.26a) 


EK±  =  pv+v (1-g ) n  1/ar(l^l/a)  =  pv+(l-g)EX1 (o)  (3  2g) 

where  EX^^°^  is  the  value  of  EX^  for  x^  =  0  given  by  Eq.  (3.17). 
Similarity 


ex2  =  gv  +  (1-p)ex2 


(o) 


(3.30) 


where  EX2  is  given  by  Eq.  (3.19).  The  expected  value  of  the  interval 
(X2~X1)  according  to  Eqs.  (3.9),  (3.29)  and  (3.30) 

E  (X2-Xl)  =  (1-P)[eX2(0)-  EX/0']  (3.31) 

while 

EX  =  pv  +  v(l-p)  T(l+l/a)  (3.32) 


Because  of  Eq.  (3.26a) 

Var  \ 

and 


Var  Y, 
k 


(3.33) 
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Var 


<xk+rV  =  var  <Yk+rV 


(3.34) 


2 

In  Fig.  2  the  ratios  EX^/EX,  E(X  -X^/EX,  Var  X^/tEX)  and 

Var  (X2~X^) / (EX) ^  are  plotted  as  functions  of  the  coefficient  of 

variation  V  of  the  distribution  for  x  =  0  and  X  =  0.1  v  (6=0.1) 
x  o  o 

and  sample  sizes  n  =  3,  20,  50,  200  and  1,000.  The  effect  of  the 

existence  of  a  minimum  life  is  quite  pronounced  although  the  assumed 

value  of  3=0.1  is  moderate.  That  a  value  of  this  order  of  magnitude 

can  be  expected  even  within  the  range  of  short  and  moderate  fatigue 

lives  is  suggested  by  Fig.  3  which  shows  values  of  x  =  N  evaluated 

10  °  ° 

from  fatigue  tests  on  aluminum  specimens  and  fitted  by  a  straight 
line,  compared  with  the  relation  Nq=  0 . lv  „  0.11N,  assuming  that  for 
3  <  a  <  4  the  relation  N  =  vr (1+1/a)  ~  0.9v  is  a  good  enough  approx¬ 
imation.  If  the  reliability  analysis  is  performed  under  the  convenient 
assumption  x^  =  0(3=0),  the  actual  existence  of  a  minimum  life  provides 
a  significant  additional  margin  of  safety.  If  the  load  spectrum  is 
specified  in  terms  of  the  percentages  p^  of  time  during  which  the  struc¬ 
ture  is  subject  to  load  amplitudes  of  certain  (constant)  intensities 
S^,  each  of  which,  if  acting  alone,  would  produce  failure  at  a  specified 

mean  time  EX.,  and  if  the  distributionsof  the  times  X.  are  assumed  ex- 
i  i 

tremal  with  characteristic  values  v.,  x  .  =  0  and  the  same  coefficient  of 

.  i  oi 

variation  (a^=  a)  >  the  reliability  function  under  the  assumption  of  no 
interaction  between  load  intensities  is  of  the  form 


P*  (x) 


p.x 


a 


a 


"i exp  (-  v~  ) -  exp  [-E  (vr)  ] 


„  (3.35) 

p .  a 

—a  r  i 

where  w  =  E( — )  :  w  is  the  characteristic  time  to  failure  under  the 

i  v. 

combined  loading.  For  a  =  1.  Eq.  (3.35)  implies  linear  damage  accumu¬ 
lation.  In  the  case  of  interaction  between  load  intensities  the  values 
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are  reduced  by  interaction  factors  uk>1  the  values  of  which  increase 
with  decreasing  load  intensity.  In  this  case  the  characteristic  value 
w  is  obtained  from  the  relation 


-a 

w 


(3.36) 


As  a  result  of  this  interaction  the  damaging  effect  of  low  load  inten¬ 
sities  is  considerably  magnified.  The  expected  time  to  first  failure 
in  a  sample  of  size  n  as  well  as  the  other  order  statistics  can  be 
obtained  from  Eqs.  (3.17)  to  (3.20)  or  from  the  diagrams  in  Fig.  2. 

Eq.  (3.7)  cannot  be  evaluated  in  closed  form  for  the  reliability 
function  Eq.  (2.17).  Numerical  evaluation  within  the  range  of  validity 
of  its  approximation  of  Eq.  (2.10)  (x^u)  <  1.0  for  the  exponents  a  = 

n  +  1  =  2,  3,  4  the  parameters  9  =  10,100  and  1,000  and  the  sample  size 
n  =  20,  50,  200  and  1,000  produces  the  values  EX^  in  terms  of  u  pre¬ 
sented  in  Table  1: 


Table  1. 


10 

100 

1000 

a 

0.362 

2 

20 

0.381 

3 

0.392 

4 

0.177 

0.836 

2 

50 

0.186 

0.820 

3 

0.190 

0.811 

4 

0.045 

0.361 

2 

200 

0.049 

0.380 

3 

0.049 

0.392 

4 

0.007 

0.088 

0.568 

2 

1000 

0.007 

0.096 

0.582 

3 

0.007 

0.099 

0.591 

4 
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The  expected  time  to  the  first  chance  failure  for  the  reliability  func¬ 
tion  Eq.  (2.8)  with  v  =  cpu  according  to  the  first  Eq.  (3.4)  is 


EX 


1C 


n  n 


(3.37) 


The  values  of  the  ratio  (EX^/EX.^)  are  presented  in  Fig.  4  as  a 
function  of  EX.^.  It  can  tie  seen  that  these  values  are  practically 
independent  of  n  and  of  a.  For  values  EX  <  0.2u  the  ratio  is  suf- 

-L 

ficiently  close  to  unity  to  justify  the  conclusion  that  under  the  as¬ 
sumptions  resulting  in  the  reliability  function  Eq.  (2.17)  fatigue  re¬ 
duces  the  expected  time  to  first  chance  failure  only  when  cp  <  0 . 2n  or 
u  <  5v/n. 

IV.  'Fatigue  Sensitivity 

The  reliability  estimate  of  a  structure  is  strongly  affected  by 
its  fatigue  sensitivity.  For  fatigue  insensitive  structures  the  re¬ 
liability  estimate  depends  on  adequate  statistical  correlation  of  the 
carrying  capacity  of  the  structure  with  a  spectrum  of  extreme  loading 
conditions,  such  as  a  tlunderstorm  or  an  extremal  maneuver  load  spectrum.  The 
statistical  varidDility  of  the  initial  carrying  capacity  has  recently 
been  established  on  the  basis  of  evaluation  of  industrial  test  data 
on  full-scale  structural  elements^  and  spectra  of  extreme  loading 
conditions  as  well  as  their  frequency  ratio  within  the  total  operation¬ 
al  conditions  can  be  reasonably  well  constructed  on  the  basis  of  extra¬ 
polation  from  flight  records  obtained  under  extremal  conditions.  The 
reliability  function  of  the  structure  for  ultimate  load  failure  can 
therefore  be  derived  on  the  basis  of  the  assumed  chance  character  of 

this  failure;  its  return  period  v  =  P  ^  is  obtained  by  evaluation  of 

F 

the  probability  of  failure  of  the  structure  from  the  correlation  of 

12 

carrying  capacity  and  load-spectrum. 

When  the  structure  is  fatigue  sensitive  the  dependence  of  the 
momentary  carrying  capacity  of  the  structure  on  its  operational 
history  destroys  the  possibility  of  this  simple  correlation.  The  ob- 
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servation  that  under  variable  .Load  intensities  fatigue  failures  in 
tests  as  well  as  in  operation  occur  at  the  highest  load  intensities 
justifies  the  assumption  that  the  chance  character  of  the 
failure  process  can  be  retained;  however ,  the  reduction  of  the  carry¬ 
ing  capacity  as  a  function  of  the  load-history  must  be  introduced. 

This  leads  to  a  reliability  function  of  the  type  of  Eq.  (2.10)  the 
actual  use  of  which  is  severely  limited  by  the  analytical  difficulties 
of  its  use.  The  proposed  approximation  Eq.  (2.17)  makes  its  use  over 
a  limited  range  possible. 

With  respect  to  this  reliability  function  the  range  of  fatigue 

insensitivity  has  been  delimited  by  the  values  of  the  expected  time 

to  first  chance  failure  0  <  EX  <  0  2u 

1C 

Thus,  for  a  fleet  size  of  n  =  50  the  condition  for  a  fatigue  in¬ 
sensitive  design  is  u  >  0 . lv  while  for  n  =  1000  u  >  0.005v.  Intro¬ 
ducing  as  an  example,  EXlc  =  2xl05  hrs.for  n  =  1000  as  an  admissible 
design  value,  the  expected  time  to  chance  failure  EX^,  =  v=  2x10®  hrs; 
a  fatigue-insensitive  design  therefore  requires  u  >  106  hrs.  For  the 
same  value  EX^^  in  a  smaller  fleet  of  n  =  50  for  which  v=  10^  hrs, 
the  requirement  for  fatigue  insensitive-design  is  again  u  >  10®  hrs. 
Assuming  therefore  that  the  design  limit  load,  considered  as  the  high¬ 
est  load  intensity  of  the  fatigue  spectrum,  is  determined  by  the  re- 
quirement  u  >  10  hrs,  independently  of  fleet  size,  the  design  ultimate 
load  intensity  should  be  the  higher  the  larger  the  fleet  size.  This 
requirement  obviously  interferes  with  the  constant  ratio  between  ulti¬ 
mate  and  limit  load  on  which  conventional  design  is  based. 

If,  for  the  sake  of  more  convenient  analytical  manipulation,  the 
reliability  function  of  a  fatigue  sensitive  structure  is  approximated 
by  the  extremal  distribution  Eq.  (2.5)  alone,  while  Eq.  (2.8)  is  the 
reliability  function  of  the  undamaged  structure  for  ultimate  load 
failure,  the  range  of  fatigue  insensitivity  is  delimited  by  the  con¬ 
dition  that  the  expected  time  to  first  ultimate  load  failure 
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Exic  be 


* 


smaller  than  the  time  to  first  failure  EX1F  of  the  fatigue  sensitive 


structure.  Considering  Eqs.  (3.17a)  and  (3.4)  the  ratio  between  EX 
and  EXlc  considered  as  inverse  fatigue  sensitivity  factor 


IF 


fl  1  =  EXiF/EXic  =  n  1/aEXp/vn  1  =  (EXF/EXc)n(1  1/a)  (4.1) 

since  v  =  EXc .  The  values  of  are  given  in  Table  2  for  n  =  3,20, 

50,  200,  1000  and  a  =  2,  3,  4. 


Table  2 


\  n 

a  \ 

3 

20 

50 

200 

1000 

2 

1.7 

4.5 

7.1 

14.2 

31.6 

3 

2.1 

7.7 

13.6 

34.2 

100 

4 

2.3 

9.4 

18.9 

53.8 

179 

These  values  are  the  ratios  between  the  design  values  (EX  /EX  )  neces- 

>  C  F 

sary  to  ensure  that  EX  =  EX  ,  a  condition  that  might  be  used  to  de- 
signate  a  structure  as  fatigue-insensitive  (f  <  1).  When  the  ratio  is 
higher  the  factor  f^  >  1  and  the  structure  becomes  fatigue-sensitive 
since  the  expected  time  to  the  first  fatigue  failure  is  shorter  than 
that  to  the  first  chance  failure.  For  n  =  1000  at  a  =  4  and  for  n  = 

50  at  a  between  2  and  3  these  ratios  are  fairly  close  to  those  ob¬ 
tained  from  Eq.  (2.17)  which,  however,  are  practically  independent  of 
a .  Nevertheless  the  fact  that  the  right  order  of  magnitude  of  these 
raitios  in  the  range  of  the  shape  parameter  2  <  a  <  4  associated  with 
the  range  of  significant  values  of  the  scatter  of  fatigue  life  0.14  < 

cr(log  X  )  <  0.19  is  obtained  from  Eq.  (4.1)  might  justify  the  use  of 
r 

this  equation  for  the  illustration  of  the  significance  of  the  concept 

of  fatigue  sensitivity  in  an  integrated  design  procedure. 

If  it  can  be  assumed  that  an  estimate  of  EX  under  a  cyclic  load 

F 

spectrum  p^(S^)  is  obtainable  from  a  modified  linear  damage  accumulation 
rule 
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I 


p.  in  EX 
i _ P_ 

^XiF 


=  k 


(4.2) 


where  EX  is  the  expected  time  to  failure  in  hours  at  the  constant 

lr 

load  amplitude  taken  either  from  an  average  S-N  diagram  by  con¬ 

verting  the  mean  number  of  cycles  N.  =  rnEX.  into  hours  with  the  aid 

1  IF 

of  the  expected  number  m  of  cycles  per  hour,  in  which  case  k  <  1,  or 
from  an  S-N  diagram  corrected  for  stress-interaction13  in  which  case 
k  =  1,  so  that 


P, 


-1 


EX_ 


=  k  f  E  — — ] 
L  i  EXiFJ 


(4.3) 


the  design  ratio 

£  =  EXc/EXf  =  [  E  ]  /  kPF  (4.4) 

i  iF 

where  Pp  is  the  probability  of  ultimate  load  failure  per  hour  of 
flight.  The  values  p^,  EX^p  and  k  are  determined  by  the  fatigue  load 
spectrum  and  the  fatigue  behavior  of  the  material  in  the  structure, 
while  combines  the  effects  of  ultimate  load  spectrum  and  ulti¬ 

mate  carrying  capacity  of  the  structure.  The  "design  ratio"  f  thus 
reflects  the  combined  principal  design  parameters.  Since  Eq.  (4.1) 
can  be  written  in  the  form 


f  =  f  n(1  -  1/a) 
1 


(4.5) 


it  takes  the  form 

f  =  n0-  -  1/a) 


(4.6) 


for  the  critical  value  =  1.  The  structure  is  fatigue  sensitive 

when  f  >  1  and  therefore  f  larger  than  the  values  tabulated  in 
Table  2  and  fatigue  insensitive  when  f  <  1  so  that  f  is  smaller 
than  these  values . 


A 
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When  a  structure  has  been  designed  for  ultimate  load  with  a 


probability  of  failure  for  a  single  load  application  p„  and  with 

r 

an  associated  probability  of  failure  per  hour  Pp  =  mpp  =  [EXC]-1, 
where  m  is  the  number  of  load  application  per  hour  in  the  range  of 
the  ultimate  load  spectrum,  Eqs.  (4.4)  and  (4.5)  can  be  used  to  deter¬ 
mine  its  fatigue  sensitivity  factor  f^  when  the  fatigue  load  spectrum 
and  the  fatigue  performance  of  the  material  in  the  structure  subject 
to  this  spectrum  are  known;  conversely,  the  maximum  intensity  of  the 
fatigue-load  spectrum  can  be  determined  for  which  the  condition  f^=l 
is  satisfied.  Thus,  for  instance  for  a  fleet  size  n  of  structures 
designed  for  ultimate  load  with  probability  pp  the  condition  of  fa¬ 
tigue-insensitivity  f i  =  1  takes  the  form  f  s  n^  or 


r 

L 

i 


Pi 


EX 


in'1’ 


(4.7) 


±F 


Assuming  an  average  of  20  gusts  per  mile  of  flight  in  thunder¬ 
storm  turbulence  all  of  which  are  of  an  intensity  that  can  be  con¬ 
sidered  to  fall  within  the  extreme  range  from  which  the  ultimate  load 
intensity  is  extrapolated,  and  assuming  further  that  flight  through 

thunderstorm  turbulence  makes  up  about  0.1  percent  of  operational 
14 

flight  time  with  average  velocity  through  turbulence  of  400  miles, 

the  expected  number  of  applications  of  loads  belonging  to  the  ultimate 

load  spectrum  is  m  =  8  gusts  per  hour  of  (total)  flight  time.  For  n  = 

1000,  a  =  4,  p  =  10“®  and  k  ~  0.5  Eq.  (4.7)  expresses  the  condition 
F 


L  EX . 

IF 

i 


< 

=  179  x  0.5  x  8  x  10 


-8 


7.15  x  10 


For  a  fatigue  load  spectrum  of  five  load  levels:  S  with  p  = 

5  4 

10-4  and  EXgF  =  Ng/m,  S4  with  P4  =  10-3  and  EX4p  =  10N5/m,  S3  with 
P3  =  10”2  and  EX3  =  100N5/m,  S2  with  p2  =  10-1  and  EX2p  =  1000N5/m 
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104N5/in  the  left  side  of  Eq.  (4.7) 
4.889  x  10  4fnN^  ^ 

The  number  m  can  be  estimated  from  the  assumption  of  an  average  of  0.1 
gusts  that  produce  fatigue  damage  per  mile  of  normal  flight  at  400  miles 
per  hour  or  m  =  40  such  gusts  per  hour.  Hence  fatigue  insensitivity  of 
the  design  requires  that  the  number  of  cycles  to  failure  at  the  highest 
load  level  S5  should  be  N5  ~  2700  cycles  or  EX5  F  =  68  hrs.  In  a  total 
estimated  life  of  25,000  hrs.  this  load  intensity  is  actually  applied 
during  2.5  hours  or  an  expected  number  of  100  times:  S,-  is  therefore 

D 

only  somewhat  lower  than  the  "limit  load"  defined  by  an  expected  fre¬ 
quency  of  occurence  of  about  10,  if  the  structure  is  to  be  fatigue  in¬ 
sensitive  with  respect  to  the  first  failure.  This  conclusion  would 

provide  some  theoretical  support  to  the  rule  of  thumb  used  in  in- 
15 

dustry  that  a  structure  or  structural  part  will  be  fatigue  insensi¬ 
tive  in  operation  if  it  can  support  about  3000  repetitions  of  the  limit 
load. 

The  fatigue  sensitivity  design  factor  f  is  a  function  of  n  accord¬ 
ing  to  Eq.  (4.5).  Therefore  the  fatigue  sensitivity  factor  ^  also 
changes  with  fleet  size.  For  the  considered  design  conditions  the 
ratio  which  is  ^  =  1  for  n  =  1000  is  increased  for  n  =  50  by  100/13.6 
=  7.3  so  that  EX^p  =  0.136EX^:  the  smaller  the  fleet  the  more  fa¬ 
tigue  sensitive  it  is  with  respect  to  first  failure.  This  only  means 
that  the  smaller  the  fleet  the  more  likely  it  is  that  the  first  failure 
is  not  an  ultimate  load  failure  but  a  fatigue  failure.  If  fatigue  sen¬ 
sitivity  were  defined  with  respect  to  the  ratio  of  expected  times  f  = 

1.0  rather  than  with  respect  to  f^,  the  fatigue  sensitivity  f-^  = 

n  (1  1/a) <  q.  the  structure  WOuld  always  be  highly  fatigue-insensitive. 


and  S1  with  p1 
becomes 


=  0.889  and  EX 


IF 
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EXiF 
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For  this  condition,  however,  the  highest  load  level  Sr  of  the  con- 

5 

sidered  spectrum  and  n  =  1000  would  be  associated  with  N  ~  5  x  105 
and  thus  be  far  below  the  limit-load. 

Definition  of  fatigue  sensitivity  by  the  ratio  of  the  risk  func¬ 
tions  associated  with  Eqs.  (2.5)  and  (2.8)^ 

a-1 


r  v  t 

=  —  =  a  —  C — ^ 

ru  VF  VV 


since 


_  jsl  r^\ 
U  / 


a-1 


and  r  =  v 
U 


-1 


(4.8) 


(4.9) 


F  F 

can  be  related  to  the  ratio  f1  with  the  aid  of  Eqs.  (3.21)  and  (3. 17) a 


a  a— 1 

t  \ 


f  =  af[r(l  +  l/a)]  (^) 


(4.10) 


Introducing  t  =  EX^  =  EXc/n 


and  therefore 


f*  =  a^[r(l  +  1/a)  j  n  (a_1) 

a 

f'  =  +  l/a)] 


(4.11) 


(4.12) 


Hence  for  f^  =  1.0 


a 

f 1  =  a[r(l  +  1/a) j  >  1 


(4.13) 


The  time  t*  at  which  the  risks  rf  and  ry  are  equal  which  has  been  pro- 
posed  as  a  safe  life  of  operation-1-  is  obtained  from  Eq.  (4.10)  with 
f  =  1 

t*  =  EXp  a_iy  (af)_1  [T  (1+1/a)  ]  _a  =  EXcf"  a/a_1  a_iy*-1  [T  (1+1/a)  ]  ~a 

(4.14) 
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Since 


-a/a-1  _  1  f  -a/a-1 

f  n  fl 


the  expected  time  to  equal  risk  of  failure 


_  -a/a-1  a-1  /  -lr_,_  .  .  .  , 

t*  =  EX1C  f±  y  a  [T  (1+1/a)] 


-a 


(4.15) 


For  a  design  with  f^  =  1  this  time  t*  <  EX^C  =  EX-^F;  for  2  <  a  <  4  the 

ratio  0.63  <  (t*/EX  )  <0.73.  Thus  the  fatigue  sensitivity  criterion 

IF 

=  1  produces  a  time  to  first  failure  which  differs  only  slightly  from 

the  time  at  which  the  failure  risks  are  equal.  This  time  t*  ~  2/3  EX^F 

can  thus  be  considered  as  a  fair  estimate  of  safe  life. 

Obviously  the  consideration  of  f^  as  the  constant  ratio  f =  EX-^/ 

EXip  is  only  an  approximation  to  the  more  rigorous  approach  in  which 

X.  and  X-,  -  are  introduced  as  random  variables  so  that  the  fatigue  sen- 
1C  -L-c 

sitivity  factor  ^  as  a  quotient  of  two  random  variables 

(4.16) 

is  itself  a  random  variable.  It  has  been  assumed,  so  far  that  the  ex¬ 
pectation  E  f ^  is  closely  enough  approximated  by  the  ratio  of  the  ex¬ 
pectations  of  Xlc  and  Xlp.  But  unless  the  distribution  of  f1  (x)  is 
known  no  probability  statement  can  be  associated  with  the  criterion 
of  fatigue-insensitivity  f  (x)  =  1. 

The  density  function  of  the  time  to  the  first  failure  pfX-j^)  is 
given  by  Eq.  (3.6)  for  any  given  distribution  of  times  to  failure  P(x). 

Thus,  for  instance,  for  the  distribution  function  Eq.  (2.5) 

„  ^  a-1  ,  ,  .a 


fl(x)  =  xlc/x1F 


P(X1F> 

=  PiF<x)  =  n 

for  a  = 

1,  becomes 

P<Xic> 

=  Plc(x)  =  n 

' 


(4.17) 


(4.18) 
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The  distribution  of  the  quotient  f^  =  (xic/xip)  is  therefore 
00 

p(fl>  =  J  Plc(xf1)plF(x)dx  = 


17 


(4.19) 


■  1  -  J na  Or)  exp  f- n  br  fi  QD  *  (t1)  ]}  d  (f-) 

F 


and 


°°  _  v  a  v  a 

\  P  2  F  /x"\  r  r  F  *  /x  ~\  ,  /"x  ,  /x  n 

P'fi>  =  J  n  —  “  OtO  “P  f-  n  L~  fi  Or; +  Cw  J]  d  0 ~J  <4-20) 


o  C  F 


The  form  of  the  distribution  and  of  the  density  function  can  only  be 
evaluated  numerically.  This  has  been  done  for  a  =  3, 4  and  n  =  3,  20 ,  50, 

200  and  1000  and  the  distribution  functions  P(f^)  are  plotted  in  Fig.  5 
on  the  scale  f^(v^/v  ).  However  the  expectation  E{f^}  and  the  variance 
a  ^  can  be  obtained  in  closed  form.  Thus  after  some  manipulation 

00 

E{fl}  =J  frP(fl)dfl  =  “  n"(1"1/a).  r  (1-1/a)  (4.21) 

o  F 

Considering  that  v  =  EX  and  v  =  EX„ [T (1+1/a) ]  ^  this  expression  can  be 

C  F  * 

written  in  the  form 

E^}  =  f  n  "  (1"1/a)r  (1-1/a)  .T  (1+1/a)  (4.22) 

of  which  Eq.  (4.5)  is  an  approximation. 

The  product  of  the  gamma  functions  in  Eq.  (4.22)  for  a  =  2,  3  and 
4  is,  respectively  (1.77x0.88)  =  1.55,  (1.35x0.89)  =  1.20  and  (1.225x0.91) 

1.10;  the  approximate  values  of  Table  2  are  the  closer  to  the  correct 
values  of  Eq.  (4.22)  the  larger  the  parameter  a.  For  n  =  1000  and  a  =  4 
used  in  the  illustration  of  the  procedure  the  correct  ratio  is  (179/110)= 
163,  which  produces  the  correct  value  of  N,-  =  3000  cycles  or  EX<_F  =75  hrs 
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The  variance 


°F1  =  E{fl2}  "  lEtV]2 


(4.23) 


Evaluating 


=  J M  „  -Ci-Vta)]2r(3, ,r(1_ 

O  -  p  J 


2/a) 


(4.24) 


valid  for  a  >  2.  The  variance 

2 


fl 


2  f  v  -(l-l/a)i  r  .  o 

=  L7"  n  { 2T (l-2/a) -  [f  (1-1/a)]2} 


=  [E{f  }]2  f  lI.Cl-2/ai,  _  . 

L  i  [r  (1-1/a)  ]  2  1i 


(4.25) 


The  coefficient  of  variation  of  is  therefore 
r  2T  (1-2/g)  . 1/2 

Vfl  MT  (1-1/a)] 2  for  a  >  2  (4.26) 


which  for  a  =  3  and  4  gives  the  values  vfl  =0.68  and  0.41.  While 
these  coefficients  are  rather  large,  they  are  associated  with  scatter 
values  of  fatigue  life  that  are  at  the  upper  limit  of  those  observed 
in  tests  of  structural  parts  and  therefore  higher  than  those  to  be 
expected  for  full-scale  structures.  Moreover  the  assumption  of  an 
extreme-value  distribution  of  fatigue  lives  for  convenient  analysis 
m  closed  form  is  quite  unfavorable  in  view  of  the  fact  that  existing 
observations  can  be  reasonably  well  fitted  by  logarithmic-normal  dis¬ 
tribution.  Thus  it  appears  that  in  reality  the  coefficient  of  varia¬ 
tion  of  might  be  smaller  than  suggested  by  the  above  values  and  that 

therefore  a  design  with  an  expected  value  0.2  <  f  <  0.4  may  produce 
structures  that  can  be  classified  as  fatigue  insensitive  with  a  suf¬ 
ficiently  high  probability.  Thus,  for  instance,  it  appears  from  Fig. 

4  that  even  under  the  assumptions  made  a  design  with  expected  E[f^]  = 
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0.3  for  a  -  4  and  would  exceed  f±  =  1  only  with  a  probability  of 
about  0.10  and  f  =  2  with  0.01. 

Obviously  in  actual  design  the  problem  arises  whether  a  fa¬ 
tigue-insensitive  design  can,  in  fact,  be  attained  with  the  ma¬ 
terial  used;  in  other  words  whether  the  operational  fatigue  load 
spectrum  applied  to  the  structure  produces  dimensions  which  are 
compatible  with  those  arising  from  the  ultimate  load  design  with 
a  prescribed  probability  of  failure  Pp.  For  a  given  material  and 
operating  conditions  such  design  may  not  be  compatible  with  con¬ 
ventionally  assumed  ratios  between  ultimate  and  limit  load;  if,  as 
a  result  a  fatigue-sensitive  design  with  f±  »  l  has  to  be  accepted 
the  necessary  provisions  for  fail-safe  construction  must  be  made 
which  may  not  be  unnecessary  in  a  fatigue-insensitive  design.  The 
estimate  of  the  safe  life  of  such  a  structure  with  a  specified  re¬ 
liability  must  be  made  on  the  basis  of  the  distribution  function 
p(f1)-  71113  subject  will  be  dealt  with  in  a  future  report. 
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Extremal  Gamma  Log  normal 


(l)d  ‘(DM 


Fig .  1  Comparison  of  frequency  and  hazard  functions  of 
Gamma,  Extremal  and  Log-Normal  distributions. 
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time, 


Fig.  2a  Ratio  EX-^/EX  versus  coefficient  of  variation  V 
for  xQ  =  0  and  xQ  =  O.lv  and  n  =  3,  20,  50,  200 
and  1000. 
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E(X2-X,) 


Fig.  2b  Ratio  E  (X2—  Xj_)EX  versus  coefficient  of  variation 
Vx  for  xQ  =  0  and  xQ  =  0 . lv  and  n  =  3,  20,  50, 
200  and  1000. 
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('x)dVA 


Fig.  2c 


Ratio  Var  X^/ (EX)  versus  coefficient 

Vx  for  xo  =  0  x  =  O.lv  for  n  =  3. 

and  1000. 


of  variation 
20,  50,  200 
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Fig.  2d  Ratio  Var (X2-X^) / (EX) ^  versus  coefficient  of  varia- 
tion  Vx  for  xQ  =  0  and  x  =  0 . lv  for  n  =  3,  20.  50. 
200  and  1000. 
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N 


Fig.  3  Minimum  life  as  function  of  mean  life  N  (Ref.  10) . 
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EX./EX 


Fig.  4  Ratio  EX1/EXlc  versus  EXlr  for  reliability  function 
Eq.  (2.17). 
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Distribution  function  of  the  fatigue  sensitivity 
factor  P(f^)  plotted  versus  fi(vp/vc)  f°r  a  =  3 
and  4  and  n  =  3,  20,  50,  200  and  1000. 
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Fig.  5 
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